Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Relativist ic Corrections to the Thermal Sunyaev-Zel''''dovich Power Spectrum

We present a quantitative estimate of the relativistic corrections to the thermal SZ power spectrum produced by the energetic electrons in massive clusters. The corrections are well within 10% for current experiments with working frequencies below v < 100 GHz, but become non-negligible at high frequencies v > 350 GHz. Moreover, the corrections appear to be slightly smaller at higher l or smaller angular scales. We conclude that there is no need to include the relativistic corrections in the theoretical study of the SZ power spectrum especially at low frequencies unless the SZ power spectrum is used for precision cosmology.     

Towards a free–free template for CMB foregrounds

A full-sky template map of the Galactic free–free foreground emission component is increasingly important for high-sensitivity cosmic microwave background (CMB) experiments. We use the recently published Hα data of both the northern and southern skies as the basis for such a template.
The first step is to correct the Hα maps for dust absorption using the 100-μm dust maps of Schlegel, Finkbeiner & Davis. We show that for a range of longitudes, the Galactic latitude distribution of absorption suggests that it is 33 per cent of the full extragalactic absorption. A reliable absorption-corrected Hα map can be produced for ∼95 per cent of the sky; the area for which a template cannot be recovered is the Galactic plane area  | b | < 5°, l = 260°–0°–160°  and some isolated dense dust clouds at intermediate latitudes.
The second step is to convert the dust-corrected Hα data into a predicted radio surface brightness. The free–free emission formula is revised to give an accurate expression (1 per cent) for the radio emission covering the frequency range 100 MHz–100 GHz and the electron temperature range 3000–20 000 K. The main uncertainty when applying this expression is the variation of electron temperature across the sky. The emission formula is verified in several extended H  ii regions using data in the range 408–2326 MHz.
A full-sky free–free template map is presented at 30 GHz; the scaling to other frequencies is given. The Haslam et al. all-sky 408-MHz map of the sky can be corrected for this free–free component, which amounts to a  ≈6  per cent correction at intermediate and high latitudes, to provide a pure synchrotron all-sky template. The implications for CMB experiments are discussed.     

三峡重庆库区深部地球物理特征与断裂构造

为了深入研究三峡重庆库区岩石圈动力学特征及其对断裂构造活动的控制和影响．并为之提供基础资料，为三峡重庆库区地震、地质灾害的监测与防治提供基础依据，在已有地球物理资料的基础上，从综合地球物理研究角度出发，通过实测地震测线资料的再解释．采用新的处理技术方法，对本区的东西向主剖面和南北向支测线的地震测深资料进行二维射线追踪处理、Pg波成像；选用场分离技术、位移数字成像技术重新处理了重力和航磁资料，通过联合反演来建立深部二维构造剖面，对剖面所揭示的基底构造特征和地壳结构特征、主要断裂构造特征以及莫霍面的起伏特征进行了精细分析和细致研究．从地球物理平面场特征出发建立了岩石圈构造三维框架。研究结果表明，这样的研究思路准确、方法得当．结论可靠；沿该剖面，把可解译的断裂分为Ⅲ级：Ⅰ级为超岩石圈断裂；Ⅱ级为壳断裂；Ⅲ级为盖层断裂。依据地球物理特征，准确揭示了库区的断裂构造特征，达到了预期效果。     

Acceleration of anomalous cosmic rays at the heliospheric termination shock

Astronomy Letters - The acceleration of anomalous cosmic rays (ACRs) at the heliospheric termination shock and their influence on the shock structure and location are analyzed in terms of a...     

Sequences in the Hardness Ratio-Peak Energy Plane of Gamma-Ray Bursts

The narrowness of the distribution of the peak energy of the vFv spectrum of gamma-ray bursts (GRBs) and the unification of GRB populations are great puzzles yet to be solved. We investigate the two puzzles based on the global spectral behaviors of different GRB populations, the long GRBs, the short GRBs, and the X-ray flashes (XRFs), in the HR - Ep plane (HR the spectral hardness ratio) with BATSE and HETE-2 observations. It is found that the long GRBs and the XRFs observed by HETE-2 seem to follow the same sequence in the HR-Ep plane, with the XRFs at the low end of this sequence. We fit the sequence by a universal Band function, and find that this sequence is mainly defined by the low energy index α, and is insensitive to the high energy index,β. With fixed β=-5, a best fit is given by α=-1.00 with X^2min/dof=2.2. The long and short GRBs observed by BATSE follow significantly different sequences in the HR-Ep plane, with most of the short GRBs having a larger hardness ratio than the long GRBs at a given Ep. For the long GRBs a best-fit yields α=-0.30 and β=-2.05. For the short GRBs, a best fit gives α=-0.60 with x^2min=1.1 (withβfixed at-2.0 because it is numerically unstable). The α value for the short GRBs is significantly greater than that for the long GRBs. These results indicate that the global spectral behaviors of the long GRB sample and the XRF sample are similar, while that of the short GRBs is different. The short GRBs seem to be a unique subclass of GRBs, and they are not the higher energy extension of the long GRBs.     

Using bright ellipticals as dark energy cosmic clocks

The age of the Universe has been increasingly constrained by different techniques, such as the observations of type Ia supernovae (SNIa) at high redshift or dating the stellar populations of globular clusters. In this paper, we present a complementary approach using the colours of the brightest elliptical galaxies in clusters over a wide redshift range  ( z ≲ 1)  . We put new and independent bounds on the dark energy equation of state parametrized by a constant pressure-to-density ratio   w _Q  and by a parameter (ξ) which determines the scaling between the matter and dark energy densities. We find that accurate estimates of the metallicities of the stellar populations in moderate and high-redshift cluster galaxies can pose stringent constraints on the parameters that describe dark energy. Our results are in good agreement with the analysis of dark energy models using SNIa data as a constraint. Accurate estimates of the metallicities of stellar populations in cluster galaxies at   z ≲ 2  will make this approach a powerful complement to studies of cosmological parameters using high-redshift SNIa.